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ABSTRACT 
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We made photometric and spectroscopic observations of the rapidly pulsating subdwarf B star EC 09582-1137 with the aim of 
determining the target's fundamental structural parameters from asteroseismology. This analysis forms part of a long-term programme 
geared towards distinguishing between different proposed formation scenarios for hot B subdwarfs on the basis of their internal 
characteristics. So far, secure asteroseismic solutions have been computed for 9 of these pulsators, and first comparisons with results 
from evolutionary calculations look promising. 
(~ i , The new data comprise ~ 30 hours of fast time-series photometry obtained with SUS12 at the NTT on La Silla, Chile, as well as 1 

^ |. hour of low-resolution spectroscopy gathered with EMMI, also mounted on the NTT. From the photometry we detected 5 independent 

I ' harmonic oscillations in the 135-170 s period range with amplitudes up to 0.5% of the mean brightness of the star. In addition, we 

O ' extracted two periodicities interpreted as components of a rotationally split multiplet that indicate a rotation period of the order of 

2-5 days. We also recovered the first harmonic of the dominant pulsation, albeit at an amplitude below the imposed 4 a detection 
t/^ . threshold. The spectroscopic observations led to the following estimates of the atmospheric parameters of EC 09582-1137: T^a - 

■ 34,806±233 K, log g = 5.80±0.04, and log N(He)/N(H) = - 1 .68±0.06. 
Using the observed oscillations as input, we searched in model parameter space for unique solutions that present a good fit to the data. 
Under the assumption that the two dominant observed periodicities correspond to radial or dipole modes, we were able to isolate a 

^ ■ well-constrained optimal model that agrees with the atmospheric parameters derived from spectroscopy. The observed oscillations are 

' identified with low-order acoustic modes with degree indices t = 0,1,2, and 4 and match the computed periods with a dispersion of 

' 0.57 %. Non-adiabatic calculations reveal all theoretical modes in the observed period range to be unstable, an important a posteriori 

, consistency check for the validity of the optimal model. The inferred structural parameters of EC 09582-1137 are T^g = 34,806 

! K (from spectroscopy), logg = 5.788±0.004, M. = 0.485+0.011 Mq, log (M,„,,/M.) = -4.39±0.10, R = 0.147±0.002 R^, and L = 

• ' 28.6±1.7 Lg. We additionally derive the absolute magnitude My = 4.44±0.05 and the distance d = 1460±66 pc. 

G^ 

o 

a^ 

■ 1. Introduction nario (e.g. lD'Cruz et al.lll996h . common envelope ejection, sta- 
J> ble Roche lobe overflow, and the merger of two he lium white 

V~l Understanding the formation of subdwarf B (sdB) stars is one dwarfs (modelled in detail bv iHan et aLHIpOl l200l . In princi- 

^ of the remaining challenges related to stellar evolution theory. It pjg^ t^e accuracy and relative importance of these evolutionary 

H ; IS generally accepted that the progenitors of these hot, compact channels can be tested by comparing the distributions of binary 

. objects (20,000 < < 40,000 K and 5.0 <\ogg < 6.2) lose ^^^j structural parameters derived from the diff^erent simulations 

too large a fraction of their envelope mass near the tip of the first jhose observed for sdB stars. While the binary properties can 

red giant branch (RGB) to ignite hydrogen-shell burning during ^e measured using well established techniques like radial ve- 

the asymptotic giant branch (AGB) phase. They instead settle on jo^ity measurements from spectroscopy, certain structural pa- 

the extreme horizontal branch (EHB), where they spend around rameters such as the mass and internal compositional layering 

10« years as helium-core burning objects surrounded by an inert ^an normally not be infeiTed using traditional methods. The for- 

hydrogen-rich envelope before joining the white dwarf cooling j^er may be derived in the case of an eclipsing binary, although 

track. However, the circumstances leading to the very fine-tuned g^^h systems are statistically rare. They can, however, be used 

mass loss necessary to produce an sdB star are still under debate. independently verify the mass estimates from asteroseismic 

Formation channels proposed feature single star and bmary evo- analysis, as done very successfully for the c ase of PG 1336-018 

lution in diff^erent flavours, including increased mass loss due to dVuckovic et al.l2007tlCharpinet et al.l2008l) . For its part, the de- 

helium-enrichment (e.g. | Sweigart| | 1997|) , the "hot flasher" see- termination of the internal layering, in particular the thickness of 

the hydrogen-rich envelope, is a pure product of asteroseismol- 

* Based on observations collected at the European Organisation for ogy. 
Astronomical Research in the Southern Hemisphere, Chile (proposal 

ID 080.D-0192). This study made extensive use of the computing facil- The discovery of rapid nonradial pulsators among subdwarf 

ities ofi-ered by the Calcul en Midi-Pyrenees (CALMIP) project and the B stars dKilkenny et al. 1997a) opened up the exciting possibil- 

Centre Informatique National de I'Enseignement Superieur (CINES), ity of using asteroseismology to probe the interiors of these ob- 

France. jects and determine the internal parameters. Often referred to 
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Fig. 1. Available H and He lines from the combined EMMl spectrum of EC 09582-1 137, overplotted with the best model atmosphere 
fit. The atmospheric parameters derived are indicated in the plot. 



as EC 14026 stars after the prototype, the rapid sdB pulsators 
exhibit luminosity variations on typical timescales of 100-200 
s with amplitudes of a few milli-magnitudes (mmags) and are 
found neai- the hot end of the EHB at 29,000 < T^ff < 36,000 
K. The oscillations have been modelled as /?(pressure)-mode in- 
stabilities driven by a classical /c-mechanism associated with an 
iron o pacity peak arising fr om the ionisation of K-shell elec- 
trons jCharpinet et alj [l996 b. It is noteworthy that models as- 
suming a uniform iron abundance in the envelope cannot excite 
pulsations unless the latter is artificially boosted to several times 
the solar value, an unrealistic proposition for these metal-poor 
stars. The problem was overcome by including radiative levita- 
tion and the resulting non-uniform iron abundance profile s in th e 
so-called "second-generation" models (Charpinet et al. Il997l) . 
These are able to explain the pulsations observed in real EC 
14026 stars to the point where the location of the instability strip 
can be accurately reproduced and the quantitative asteroseismo- 
logical interpretation of pe riod spect ra for individual stars has 
become possible (see Font aine et al.l i2p08; 0stensenl |2009| for 
recent reviews on the subject). To date, full asteroseismic anal- 
yses have been carried ou t for 9 rap i dly pulsating sdB stars: 
PG_0014+067_jBrassard etl] llOOlt ICharpinet et all l2005bl 
Brassard & Fontaine! l2008irPG 1047 -h003 dCharpinet etall 



20031). PG 12 19-1-534 (ICharpine t et all l2005ch . Feige 48 
(ICharpinet etall I20q5al: Ivan Grooteletal] l2008al). EC 20117- 
4014 (Randall et al.' '2006^ ), PG 1 325 -H 101 (Charpinet et al.' 
L2OO6), PG 0911+456 (RandalLeLalJ!2007), Balloon 090100001 
dvan Grootel et alJ l2008bl) and PG 1336-018 (ICharpinet etld] 
I2OO8I) . In this Paper we present the tenth asteroseismological 
analysis of an EC 14026 star. 

The target of the present study, EC 09582-1137, was 
identified in the Edin burgh-Cape (EC) survey zone 1 paper 
dKilkennv et al.lll997bh as a subdwarf B star with no spec- 
troscopic signature of a cooler companion. Infrared mea- 



surements from the 2MASS database (see the web page 
www.ipac.caltech.edu/2mass) tend to confirm this, although 
their accuracy is rather poor because of the relative faintness 
of the star (y=15.26 from the EC survey). EC 09582-1137 
was discovered to b e pulsating only relatively recently by 
iKilkennv et alj ( l2006l hereafter K2006), who uncovered two pe- 
riodicities at 136.0 and 151.2 s with amplitudes of 8.3 and 7.5 
mmags respectively. No further oscillations were detected down 
to a threshold of 1 mmag. 

In what follows, we present new spectroscopic and photo- 
metric observations of EC 09582-1137. These are used to de- 
rive atmospheric parameters and to detect oscillation frequencies 
beyond those already known (Section 2). We then provide de- 
tails of the asteroseismological analysis conducted on the basis 
of our observations, and discuss the internal parameters inferred 
(Section 3). Finally, we summarise our results and put them into 
the context of current research in Section 4. 

2. Observations & analysis 

2.1. Spectroscopy and atmospheric analysis 

We were allocated one hour of service mode observ- 
ing time with EMMI at the 3.5-m NTT located at La 
Silla, Chile (for details on the instrument see the web 
page http://www.eso.org/sci/facilities/lasilla/instruments/emmi). 
Three low-resolution spectra, each with an integration time of 
1000 s, were obtained on 24 December 2007 using Grating 
4. Each spectrum was bias and flat field corrected, extracted, 
sky subtracted and wavelength calibrated using standard IRAF 
routines. The combined spectrum has a wavelength resolution 
A/1 ~ 5 A and signal-to-noise of S/N ~ 160/pixel at the cen- 
tral wavelength Ac- 4355 A, and covers the wavelength range 
from 3413 to 5278 A. It was analysed using a detailed grid of 



S.K. Randall et al.: Asteroseismology of EC 09582-1137 



3 



non-LTE model atmospheres and synthetic spectra designed es- 
pecially for subdwarf B stars. These were computed using the 
pubhc codes TLUSTY and SYNSPEC (Hubenv & Lanz 1991 
lLanz&Hubenvlll995h . and include helium but no metals. Fig. 
[U shows the best simultaneous model fit to the available Balmer 
and helium lines, and indicates the atmospheric parameters in- 
ferred. The dip near Hy is not understood and interpreted as a 
glitch in the data, of unknown origin. As typical in sdB stars 
with such parameters, the Hell 4686 line is observed stronger 
than predict ed, which is likely due to the lack of metals in our 
models (see iHeber et"aLll2000l) . Note that for sdB stars around 
Tpff ~ 35,000 K this does not severely affect the atmospheric pa- 

PE — 



rameters deri yed (iGeier et a 



e.g. Fig. 1 of iFontaine et al 



2007 ). Comparing the latter with 



(I2008h reveals EC 09582-1137 to 



be a typical EC 14026 star in terms of atmospheric parameters, 
located towards the hot end of the instability strip. 

2.2. Time-series pliotometry and frequency analysis 

For the photometry part of the observations we were allocated 
5 nights with SUSI2 (SUperb Seeing Imager), mounted on the 
NTT at La Silla. The observing log is given in Table [T] The full 
dataset comprises ~ 30 hours of time-series photometry spread 
over 103 hours, which yields a frequency resolution of 2.7 yuHz. 
Given the relative faintness of the star, we chose to keep the filter 
slot empty with the aim of collecting as many photons as possi- 
ble. In order to reduce the readout time we employed the SUSI2 
'windowing' mode, which reads out only a pre-defined part of 
the chip, and moreover used 3x3 binning. Thanks to the small 
pixel scale of the SUSI2 imager (0.08"/pixel without binning), 
the FWHM of the stars in the images was always above ~ 3 pix- 
els, even during periods of good seeing (the seeing on the images 
ranges from 0.7"to 4", with an average of 1.3")- The resulting 
dead time between exposures is 13 s, which together with the 
exposure time of 10 s gives a cycle time of 23 s. 



Table 1. Photometry obtained for EC 09582-1 137 (2008) 



Date (UT) 


Stai-t time (UT) 


Length (h) 


10 Mar 


00:52 


00:51 


10 Mar 


06:27 


01:23 


11 Mar 


03:55 


04:20 


12 Mar 


00:40 


07:11 


13 Mar 


00:17 


07:44 


14 Mar 


00:19 


07:41 
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Fig. 2. All light curves obtained for EC 09582-1 137 with SUSI2. 
The data have been shifted arbitrarily along the x and y axes for 
visualisation purposes. From top to bottom the curves refer to the 
nights of 10, 11, 12, 13 and 14 March 2008. Details are given in 
TablelU 
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Fig. 3. Expanded view of a portion of the light curve obtained 
on 13 March 2008 in units of fractional brightness intensity and 
seconds. 



The data were reduced using standard IRAF aperture pho- 
tometry routines, except that the photometric aperture size was 
adjusted to 2.75 times the FWHM in each image. Differential 
photometry was obtained on the basis of 3 suitable comparison 
stars of similar brightness to the target that were specifically in- 
cluded in the windowed area of the CCD chip. The resulting 
light curves for all 5 nights are displayed in Figure |2] where the 
large gaps in the top two light curves are due to the telescope 
being closed because of high humidity. It is apparent that the 
longest two light curves corresponding to the last two nights are 
degraded in quality towards the end; this is a result of the high 
airmass {X ~ 2.3) at which the setting target was observed. An 
expanded version of one of the light curves is shown in Figure|3] 
where the pulsation and beating between individual frequencies 
can be better appreciated. 

We computed the Fourier transforms (FT) for each of the 
light curves individually as well as for different combinations 



of data sets. The individual nightly Fourier transforms are plot- 
ted in Figure |4] In the end, the lowest noise level was achieved 
by combining the entire data set, which is not all that surprising 
considering the relative homogeneity of the data quality. The re- 
sulting Fourier transform is displayed in the top panel of Figure 
|5]in the zoomed-in 5-9 mHz range, the remaining power spec- 
trum out to the Nyquist frequency of 21.7 mHz being consistent 
with noise. The lower curves refer to successively pre-whitened 
data, as indicated. During the process of pre-whitening, the dom- 
inant peak is identified from the FT and used as input for a least- 
squares (LS) fit to the light curve in which the amplitude and 
phase are determined. The resulting sinusoid is then subtracted 
from the light curve, and the sequence is repeated until all pe- 
riodicities down to a specified amplitude threshold have been 
extracted. The subtraction of periodicities in time rather than fre- 
quency space ensures that not only the actual FT peak, but also 
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Fig. 4. Fourier transforms for the individual light curves for each 
night in the 3-11 mHz range. The curves have been shifted ar- 
bitrarily along the y-axis for visualisation purposes, however the 
amplitude scale is the same as for the original Fourier transform. 
The locations in frequency space of the 5 independent harmonic 
oscillations extracted down to a threshold of 4 cr are marked by 
the dotted vertical lines. 

the sidelobes caused by the daily gaps in the data are removed in 
each successive FT. 

We detected a total of 6 frequencies down to an imposed 
threshold of 4 times the average noise level. In Table |2] we list 
these frequencies together with the corresponding periods, am- 
plitudes and phases, as well as the S/N derived from the ratio 
of the amplitude to the Icr noise level (measured to average 
0.0071 % over the 0-16 mHz frequency range). The periodic- 
ities are numbered according to amplitude rank, and the corre- 
sponding frequencies and amplitudes reported by K2006 are in- 
dicated where applicable. We also include two lower amplitude 
periodicities (marked by asterisks in the S/N column): the prob- 
able multiplet component /3-1- (see below), and the harmonic of 
the dominant frequency fiH. While the harmonic has a S/N of 
only 2.3, and as such falls below the threshold for credible peaks 
by any standards, it is found at exactly twice the frequency of 
the dominant oscillation, and is thus likely to be real. It has no 
part to play in the asteroseismological analysis detailed below, 
however it is interesting to note that we detect it for this sdB star 
despite the lack of significant convection in the envelope, which 
is sometimes invoked to explain the pres ence of ha rmonics in the 
pulsation spectra of other types of stars (IWull200j]) . 

The errors on the period (frequency), amplitude, and 
phase were calculated u si ng th e recipes described in 
iMontgomerv & Odonoghud d 19991 hereafter MO). These 
give uncertainties on the amplitudes and phases virtually 



identical to those obtained from the LS fit to the light curve, 
at least when applied to the full data set. When considering 
one night at a time, the more conservative MO recipes yield 
uncertainties around 10 % larger, and it is these estimates that 
we employed for the modelling of the amplitude and phase 
variations discussed below. In addition to the 8 frequencies 
listed in Table |2l two further low amplitude periodicities were 
found just above the less stringent threshold of 3cr at 135.750 
s (A ~ 0.022 %) and 132.826 s (A ~ 0.022 %). These are not 
considered reliable enough as input for the initial asteroseismic 
analysis, but are nevertheless interesting to note. 

Examining the periodicities listed in Table|2l we find that the 
period spectrum of EC 09582- 1 1 37 is dominated by two peaks at 
~136 s and ~151 s respectively. These have amplitudes an order 
of magnitude higher than the other oscillations, and correspond 
to the two frequencies detected by K2006. However, the ampli- 
tudes of both, and in particular the secondary peak are much 
higher in the earlier data. This could be partly due to the slightly 
different QE curves of the SUSI2 CCD and the UCT-CCD used 
by K2006, however upon comparison these were found to be 
quite similar. In any case, the dramatic relative change in power 
between the two peaks cannot be explained in terms of different 
CCD responsitivities. It is more likely that the amplitudes have 
changed either intrinsically or as a result of the beating between 
unresolved modes (our data have a longer baseline). The weaker 
oscillations fall below the K2006 detection threshold of 0.09 %, 
so it is impossible to say whether their amplitudes have changed 
or not. Looking at the frequencies of the two main peaks, these 
seem to be systematically lower by ~ 0.5 juHz in our data. It is 
unclear whether this is a coincidence as the potential frequency 
shift is similar to the accuracy of the measurements, or whether 
it could be a tentative manifestation of radial velocity variations 
caused by the presence of an unseen companion. 

An interesting feature of the period spectrum detailed in 
Table |2] is the apparent multiplet structure around /3. The two 
higher frequency components and /3-1- have comparable am- 
plitudes well above the 4cr limit, while /s- falls below the 
threshold. However, we nevertheless believe that it corresponds 
to a sinusoidal variation in the data rather than a spurious peak. 
In order to convince the reader of this. Fig. |6] shows the light 
curve prewhitened by the 7 other oscillations of Table |2] and 
folded onto the 143.247 s component. The individual data points 
have been combined into 8 phase bins for clarity, and the sinu- 
soid derived from the least-squares fit to the light curve has been 
overplotted. It is clear that, within the estimated errors, the points 
provide a good fit to the sinusoid, as they should if the peak is 
real. 

The most commonly invoked explanation for closely spaced 
frequency multiplets in subdwarf B stars is rotational splitting. 
It is well known that in a rotating star spherical symmetry is 
broken and the m-fold degeneracy of a mode characterised by 
degree I and radial order k is lifted. Treating rotation as a first- 
order perturbation under the assumption of solid-body rotation 
results in a frequency splitting between adjacent values of the 
azimuthal index m (where m - (-( + 1), {{ - 1), £) of 

A/^i^ (1) 

"rot 

where Proi is the rotation period and Ca is the dimensionless 
first-order rotation coefficient. The latter is normally negligible 
compared to unity for the /?-modes observed in EC 14026 stars 
(see Table |4]for the example of the optimal model). 

In the case of EC 09582-1 137, the multiplet is not evenly 
spaced in frequency, with A/3,/3- = 5.12 + 0.35 /iHz and 
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Rank 


Period 


Frequency 


Amplitude 


Phase 


S/N 


fK2Q06 






(s) 


(mHz) 


(%) 


(s) 




(mHz) 


(%) 


.h 


169.1109±0.0054 


5.91328+0.00019 


0.0445+0.0056 


49.02+3.38 


6.3 






fs 


162.4811+0.0060 


6.15456+0.00023 


0.0371+0.0056 


107.41+3.90 


5.2 






fi 


151.2420±0.0008 


6.61192+0.00004 


0.2297+0.0056 


100.26+0.59 


32.4 


6.61254 


0.69 


h- 


143.2468±0.0067 


6.98096+0.00033 


0.0258+0.0056 


10.46+4.95 


3.6* 






h 


143.1418+0.0024 


6.98608+0.00012 


0.0726+0.0056 


14.99+1.75 


10.2 






h+ 


143.0664+0.0026 


6.98976+0.00013 


0.0663+0.0056 


63.31 + 1.92 


9.3 






fi 


135.9967+0.0003 


7.35312+0.00002 


0.4739+0.0056 


14.67+0.26 


66.7 


7.35342 


0.76 


m 


67.9970+0.0024 


14.70653+0.00052 


0.0163+0.0056 


19.29+3.71 


2.3* 







* below 4(T 



A/3,/3-f = 3.68 + 0.18 yuHz. We believe that this is related 
to the multiplet not being completely resolved, even in the full 
data set. Indeed, the frequency resolution of the combined pho- 
tometry is 2.7 yuHz, of the same order of magnitude as the fre- 
quency splitting measured, and when analysing data subsets of 
2-4 continuous nights we found A/ to decrease as the time 
baseline increased. Going back to the Fourier transforms for the 
individual nights (see Fig. |4]l, we find that the 143 s structure 
manifests itself as one broad peak that appears to significantly 
change its amplitude from one night to the next. To investigate 
this in more detail, we show the nightly variation in amplitude 
and phase of the (unresolved) 143 s peak in Fig. |7] The over- 
plotted dashed (continuous) curves refer to the amplitude and 
phase variations expected from the beating action of the triplet 
structure /3-,/3,/3+ (doublet structure f-i,f-i+). While both the 
doublet and triplet scenario can account for the amplitude varia- 
tions observed relatively well, the triplet model is slightly more 
successful at matching the observed phase variations. The agree- 
ment in amplitude is by no means perfect, particularly for the 
first two nights of observation, where the data sets were some- 
what shorter than for the last three nights. We interpret this in 
terms of an insufficient resolution of the dataset: when extract- 
ing the nightly amplitudes and phases of the 143 s structure, we 
kept the period fixed at the /3 value listed in Table |2] If that pe- 
riod is poorly measured because of insufficient resolution, this 
will induce small errors in the nightly amplitude and phase mea- 
surements. Nevertheless, we believe that, together with Fig. |6] 
the results of Fig.|7]point towards the lower amplitude f-^ - com- 
ponent being real, and the 143 s complex containing at least a 
triplet of peaks. In this case one could infer an approximate ro- 
tational period of f ~2.63 days based on the mean frequency 
spacing of A/=4.4 juHz. However, given the limited resolution 
of our dataset it is of course quite possible that the structure is 
in fact an unresolved quintuplet. In this case we may be see- 
ing the m--2, m-Q, and m-+2 components of an (-2 mode, the 
frequency spacing would correspond to twice the splitting from 
rotation, and the rotation period of the star would be 5.26 days. 

It is unfortunate that we do not detect rotational splitting in 
any of the other peaks, as this would allow us to determine the 
rotational period with more certainty. While it is quite likely that 
one of the dominant periodicities is a radial mode and as such 
would not be subject to rotational splitting, at least some of the 
other modes must have higher degree indices. For the lower am- 
plitude modes fn and /s the multiplet components probably fall 
far below our detection threshold, but for f\ and fi they could be 
measurable, assuming of course that these are not radial modes. 
In this context, the 3cr oscillation at 135.750 s becomes interest- 
ing as a potential multiplet to the dominant 135.997 s periodicity. 
However, the frequency splitting is A/=13.5 yuHz, much larger 
than that attributed to rotational splitting from the /3 multiplet, so 



we believe it to correspond to an independent harmonic oscilla- 
tion with no direct connection to the strongest peak. Considering 
the case of /2= 15 1.242 s, we indeed find a very closely spaced 
component just below the 3cr level at 151.206 s, which would 
imply a frequency separation of A/=1.5 yuHz. Since the latter 
is below the formal resolution of the dataset it cannot be taken 
at face value, but it does hint at the presence of some multiplet 
structure, quite possibly with a spacing around half that found 
for /3. This would then imply { >2 for /3, however a longer data 
set is needed for confirmation. 

In summary, we detected 5 independent harmonic oscilla- 
tions as well as one peak presumed to constitute a rotationally 
split multiplet component above the Act threshold in our data. 
Imposing a lower threshold of 3 cr we additionally find a sec- 
ond split multiplet component and two lower amplitude frequen- 
cies. While the multiplet component forms part of a triplet and 
is thought to be real, the weaker independent oscillations are 
thought of as insecure detections. Finally, we uncovered the first 
harmonic of the dominant oscillation at a very low amplitude. 

3. Asteroseismic analysis 

3.1. Methodology 

For the asteroseismological analysis of EC 09582-1137 we fol- 
low the well-known f orward method described in detail by 
ICharpinet et al.l (l2005cl) . This approach is based on the require- 
ment of global optimisation, implying that all the periods ob- 
served are simultaneously matched to those computed from sdB 
star models. For the latter we employ ou r so-called "second- 
generation" (2G) models (see lCharpinet et a l. 1996, 2001), static 
structures composed of a hard ball nucleus surrounded by a more 
accurately modelled envelope extending down to a logarithmic 
depth oflog^ = log(l - M{r)/Mt) - -0.05, sufiicient for an ac- 
curate computation of the shallow /?-modes observed in the EC 
14026 stars. An important feature of the 2G models is that they 
incorporate microscopic diffusion under the assumption of an 
equilibrium having been established between gravitational set- 
tling and radiative levitation. For the case of the sdB star models 
it is iron that is assumed to be levitating in a pure hydrogen back- 
ground. Since it is the local overabundance of iron in the driving 
region that creates the opacity bump necessary for the excita- 
tion of pulsations through the /c-mechanism, including diffusion 
processes is vital when it comes to predicting whether a mode 
is unstable or not. Moreover, microscopic diffusion changes the 
stellar structure sufficiently to affect the frequencies of pulsa- 
tions and must therefore be taken into account when attempt- 
ing quantitative analyses (Fontaine et al. 2006a). The 2G models 
are specified by four fundamental input parameters: the effective 
temperature Tgff , the surface gravity log g, the total stellar mass 



S.K. Randall et al.: Asteroseismology of EC 09582-1137 



6 

Mt, and the fractional thickness of the hydrogen-rich envelope 
log^(H) = logM(H)/M». The latter parameter is intimately re- 
lated to the more commonly used quantity Mf„i,, which corre- 
sponds to the total mass of the H-rich envelopfl 

The 2G models are the input for the computation of adi- 
abatic and non-adiabatic oscillation mod es using two efficient 
codes based on finite element techniques dBrassard et alJ[T992t 
iFont aine & Brassard 1994). The first solves the four adiabatic 
oscillation equations and constitutes an intermediate but neces- 
sary step to derive estimates for the periods and eigenfunctions 
that are then used as first guesses for the solution of the full, 
non-adiabatic set of oscillation equations. The latter computes a 
number of quantities for each mode, including some that can be 
directly compared with observed quantities, such as the period 
Pfhi- 27r/crR, where o-r is the real part of the complex eigenfre- 
quency), and the stability coefficient crj (the imaginary part of 
the eigenfrequency). If ctj is positive, the mode is stable, while 
if it is negative the mode is excited, and may therefore be ob- 
servable if its amplitude is large enough. Although the aster- 
oseismological analysis is based on the more accurately com- 
puted adiabatic modes, the non-adiabatic approach can be used 
as a consistency check in the sense that the periods observed 
should also be predicted to be unstable. As is standard, the pulsa- 
tions are computed assuming perfect spherical symmetry, which 
is well justified for slowly rotating stars such as EC 09582-1137. 
This implies that each theoretical mode is (2^-Hl)-fold degener- 
ate in eigenfrequency, and all multiplet components seen in the 
observed spectrum (such as the three /? components identified 
in the present case) must be considered as a single harmonic os- 
cillation for the purpose of the asteroseismological analysis. In 
this context, a pulsation mode is then completely defined by the 
degree index i and the radial order k. 

Our approach to asteroseismology relies on a double opti- 
misation procedure that first determines and quantifies the best 
match between the set of observed periodicities and those cal- 
culated for a given model, imposing certain restrictions on the 
mode identification if necessary (see Section 3.2). Subsequently, 
we identify the model (or family of models) that can reproduce 
the observed periods most accurately within pre-defined limits 
of the 4-dimensional (Teff, log^, M*, log ^(H)) model parameter 
space. The latter is known as the "optimal" model, and corre- 
sponds to the absolute minimum of the goodness-of-fit merit 
function found by a dedicated optimisation code based on 
a hybrid Genetic Algorithm (GA, see Ch arpinet et al..,2008. for 
details). The merit function is given by 

S'-tj^(PL-P',Hf (2) 

!=l 

where is one of the n periodicities observed, P'^j^ is the 
ffieoretical period that matches it best, and cr is the global weight, 
defined as the inverse of the mode density of each model (i.e. the 
ratio of the width of the considered period window to the num- 
ber of modes in that window). Including the global weight in 
the computation of removes, at least partially, possible large 
scale biases towards models with higher mode densities, which 

' Note that the parameter M„„, commonly used in extreme horizon- 
tal branch stellar evolution includes the mass of hydrogen contained in 
the thin He/H transition zone, whereas the parameter M{H) used in our 
envelope models does not. They can be related with log [Mf„y/M,] = 
log [M(/f)/M,] + C, where C is a small positive term slightly depen- 
dent on the model parameters that can be computed from the converged 
model using the mass of hydrogen present in the transition zone itself. 



are usually found at ffie low surface gravity end of the parame- 
ter space explored. It is largely a technical choice, designed to 
prevent the GA code from allocating too many computational 
resources to a region of parameter space with unrealistic models 
that are in conflict with spectroscopic estimates of \ogg. 

3.2. Search for the optimal model 

We searched parameter space for an optimal model based on the 
5 independent harmonic oscillations listed in Table|2]as /i to /s. 
For the reasons explained above, the rotationally split compo- 
nents for fi were not included, and the lower amplitude oscilla- 
tions identified were deemed too insecure for the asteroseismo- 
logical search. As has become standard in recent asteroseismo- 
logical studies, the effective temperature was kept fixed at the 
spectroscopically determined value of Tetf =34,806 K, while the 
other parameters were confined by limits designed to sandwich 
EC 09582- 11 37 in 3-dimensional log g-M, -log q{li) parameter 
space. This yields the most accurate results because the p-mode 
frequencies in sdB stars are only weakly dependent on the effec- 
tive temperature, and the latter can generally be determined more 
accurately from spectroscopy than asteroseismology. Moreover, 
it is well known that the period spectra computed for different 
models are subject to a degeneracy in Teif - M, space, implying 
that one of the two parameters mu st be constrained in or der to 
be able to infer the other (see e.g. ICharpinet et al]|2005c ). The 
boundaries set for the surface gravity cover a generous range 
centred on the spectroscopic estimate, with 5.70 < log^ < 5.90. 
For the remaining two parameters we relied on stellar evolu- 
tion theory and th e possible ranges computed for various evolu- 
tionary scenarios (iHan et al.ll2002l I20b3i) . setting limits of 0.30 
< MJMq < 0.7 and -5.20 < log^(H) < -2.00. We computed 
modes from 50-350 s, which amply covers the observed period 
spectrum. 

It is important to mention that, although we do not have ac- 
cess to a priori mode identification from e.g. multi-colour pho- 
tometry or time-series spectroscopy, we place some constraints 
on the degree index of the observed modes in what follows. This 
is based entirely on mode visibility arguments, and assumes a 
similar intrinsic am plitude for all t he mo des considered. As il- 
lustrated in Fig. 9 of'R andall et al.l (12007'). mode visibility when 
integrated over the visible disk of the star generally decreases as 
the degree index increases for the vast majority of inclinations 
of the pulsational axis compared to the line of sight, radial and 
dipole modes dominating the amplitude hierarchy for an incli- 
nation angle / < 70°. It has moreover been shown that because 
of their specific surface geometry, modes with {=3 have an ex- 



Table 3. Models that can account well for the five secure har- 
monic oscillations detected for EC 09582-1 137. 



Rank 


Period (s) 


C k 


{ k 


/4 


169.10 


2 2 


2 1 


/5 


162.48 


4 2 


4 1 


fl 


151.24 


2 


1 2 


h 


143.14 


2 3 


2 2 


/l 


136.00 


3 


2 






Model I 


Model II 






0.146 


0.188 




AP/P 


0.41 % 


0.57 % 




log 5 


5.749 


5.788 






0.70 


0.49 




log^(H) 


-3.34 


-4.48 
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37000 36000 35000 34000 33000 32000 31000 30000 29000 -5.20 -4.80 -4.40 -4.00 -3.60 -3.20 -2.80 -2.40 -2.00 



Teff (K) log q(H) 

Fig. 8. Left panel: Slice of the S ^ function (in logarithmic units) along the Tetf - log g plane with M* and log ^(H) set to their optimal 
values of M»=0.49 and logq'(H)=-4.48. The solid rectangle represents our spectroscopic estimate of the atmospheric parameters 
for EC 09582-1137. Also indicated are the 1,2 and 3cr limits on the minimum (white dashed contours). The errors on the 
asteroseismologically derived parameters are estimated from the semi-axes of the Icr contour. Right panel: The same as the left 
panel, except that the slice of the S ^ function lies in the M» - log ^(H) plane and the other two parameters are set to their optimal 
values of refF=34,806 K and logg=5.788. 



tremely low visibility in the optical domain due to cancellation 
effects (Randall et al. 2005). Therefore, we consider only modes 
with ( - 0,1,2,4 in our asteroseismological analysis. In an ini- 
tial search, the identification of a given observed periodicity was 
left open within these constraints, however this produced a large 
number of models that were able to account for the frequencies 
observed quite well. Given that the two highest amplitude oscil- 
lations /i and /2 completely dominate the observed period spec- 
trum, we thought it justified to limit them to radial or dipole 
modes, while leaving the mode identification of the lower am- 
plitude modes open to ^ = 0,1,2,4. 

In the search domain specified, this approach identified two 
families of models able to reproduce the observed periodocities 
within the constraints in terms of /?-modefl These are listed in 
Table [3] together with the mode identification inferred. It is ap- 
parent that both models are consistent with the spectroscopic 
estimate of \ogg within the errors (Teff was fixed to the spec- 
troscopic value of 34,806 K), although Model II lies closer to 
the value measured. The two models feature a plausible mode 
identification, the lower amplitude modes being associated with 
modes of relatively high degree index. In terms of the quality of 
the period match Model I fares slightly better than Model II, but 
the fit is acceptable in both cases. 

Looking at the model parameters more closely, we find that 
Model I has a rather high mass, at the edge of the search lim- 
its defined. This alone does not exclude it as a viable model, 
since masses significantly higher and lower than the canonical 
value are possible according to binary evolution calculations 

^ A third model was identified, however it contained mixed modes, 
i.e. modes that show characteristics of both p- and g-modes. Since 
these propagate into the deeper layers of the star, they cannot be ac- 
curately modelled with the 2G models. Therefore, this model could be 
excluded. 



(iHan et al.l l2002l 120031; IYu & Lil 12009 '). However, the relative 
values of M, and log q'(H) derived are implausible from an evo- 
lutionary point of view. Regardless of how it is formed, an sdB 
star is expected to start its helium burning lifetime on a zero-age 
EHB (ZAEHB) specified by its total mass and initial metallicity. 
Its exact location on that ZAEHB is determined by the thickness 
of its hydrogen rich envelope. So at f=0, for a given metallicity, 
a precise correlation exists between the total mass, the thick- 
nes s of the H shell, and the effective temperature (for details 
see iFontaine et al.ll2006bl) . This correlation is shown in Fig. |9] 
for total stellar masses of M^IMq-QAQ, 0.45 and 0.50 assuming 
a metallicity of Z=0.003. While the high mass of M/M,=0.70 
predicted for Model I is not shown, it is clear from an extrapola- 
tion of the plot that, given the spectroscopic temperature Tgff ~ 
35,000 K one would expect a much more massive envelope than 
that indicated in Table[31 with logM(H)/M, > -2. 

The relations displayed in Fig. |9] are assumed to be valid 
to first order during most of the He-burning lifetime of sdB 
stars, but break down for individual objects that have already 
evolved away from the Z AEHB (see the unusual cas e of Balloon 
090100001 presented by Ivan Grootel et al.l l2008bl) . The loca- 
tions on the plot of the stellar parameters derived for the tar- 
gets so far submitted to asteroseismology are in accordance 
with the theoretical relations, nicely illustrating the strong in- 
ternal consistency that exists between the derived parameters. 
The one exception (marked by the outlying cross) is for Balloon 
090100001, a supposedly more evolved object. Given that EC 
09582-1137 shows both atmospheric and pulsational character- 
istics typical of EC 14026 pulsators, it makes sense to assume 
that the ZAEHB relations hold true, allowing us to exclude 
Model I as physically implausible. 

We thus retain Model II as the optimal model, and show the 
behaviour of the merit function in its vicinity in parameter space 
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Table 4. Period fit for the optimal model. 
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(s) 


(rad/s) 


(erg) 




(%) 
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103.0572 


-5.860 X 10"^ 
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3 




113.0118 


-2.057 X 10"^ 


40.575 













2 


136.0000 


134.7030 


-8.619 X 10"^ 


40.897 




+0.9537 


+ 1.2970 


111 





1 




156.8309 


-2.918 X 10"' 


41.986 


















172.1053 


-1.004 X 10"^ 


41.956 










1 


5 




101.1261 


-5.016 X 10"^ 


40.053 


0.01330 








1 


4 




110.5917 


-3.274 X 10"^ 
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0.01420 








) 


3 
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2 


151.2400 


151.1447 


-5.241 X 10"' 
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+0.0630 


+0.0953 


tJ2 


1 


1 




171.2927 
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2 


4 
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0.01655 








2 


3 
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-8.736 X 10"^ 
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2 


2 


143.1400 
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+0.6667 
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2 


1 
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169.9215 
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-0.4858 
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3 
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3 


1 
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4 


1 




271.3005 


+6.164 X 10"'2 


46.522 


-0.03263 









in Fig. [8] It is evident that the atmospheric parameters inferred 
agree perfectly with the spectroscopic estimates for \ogg and 
Tetf . As predicted, the solution is far more sensitive to the surface 
gravity than the effective temperature because of the stronger pe- 
riod dependence. In M» - log ^(H) space, the solution is well de- 
fined for both parameters, as is to be expected from their strong 
signature on the frequency spectrum. Examining the mode iden- 
tification, we see that the observed amplitude of a periodicity 
decreases with increasing degree index i, which makes the so- 
lution plausible from a mode visibility point of view. Of course, 
the mode identification was partly fixed from the outset, as f\ 
and fi were required to be radial or dipole modes. However, it 
was not obvious that the lower amplitude modes would be as- 
sociated with the less visible {-2 and {-A modes. This is par- 
ticularly interesting for the case of /s, which from the available 
data shows a triplet structure attributed to rotational splitting. If 
it is indeed a quadrupole mode as suggested by our mode iden- 
tification, we should see a quintuplet given a long enough time 
baseline. With the current data set, we are presumably able to de- 
tect only the outer m - +2 components, and cannot resolve the 
m - +\ peaks. This could explain why we do not observe rota- 
tional splitting for the higher amplitude /2 periodicity, while the 
multiplet components of the weak fn and fs peaks are probably 
lost in the noise. A longer data set would shed light on this, and 
could be exploited to confirm the mode identification and thus 
the validity of our optimal model. 



3.3. Period fit and structural parameters for EC 09582-1137 

The optimal model can account for the 5 observed periods to 
within 0.57 % (corresponding to an absolute period dispersion 
of AP=0.833 s), which is quite typical for the asteroseismologi- 
cal analyses carried out to date for EC 14026 stars. Details on the 
optimal model period spectrum can be found in Table ID where 
we list the periods computed together with their degree index { 
and radial order k, kinetic energy log E, rotation coefficient Cki 
and stability coefficient ctj. When the latter is negative, the oscil- 
lation is predicted to be excited from non-adiabatic theory. The 
fact that this is the case for all modes assigned to an observed pe- 
riodicity constitues an important a posteriori consistency check 
for the validity of our optimal model. 

The periods observed for EC 09582- 11 37 are shown next to 
their theoretical counterparts, and the absolute and relative dis- 
persion is indicated for each. A posteriori it is interesting to note 
that the two lower amplitude periodicities detected at 132.83 s 
and 135.72 s (see Section 2.2) can be respectively matched to 
the {=1, k=3 and the {=4, k-2 mode without a significant de- 
terioration of the merit function occurring. In fact, if we repeat 
our asteroseismic search including these extra two frequencies, 
we find an absolute minimum in the merit function at values that 
are compatible with the adopted optimal model within the er- 
rors. The absolute period dispersion then increases slightly, but 
the relative period fit remains unchanged. 

In order to better illustrate the quality of the period fit, we 
include a graphical representation in Fig. [TO] As has been noted 
in previous studies, the dispersion of the period fit is an order of 
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Table 5. Inferred structural parameters for EC 09582- 11 37 



Quantity 


Estimated Value 




Teff (K)T 


34806 = 


t 230 


(0.66%) 




5.788 = 


t 0.004 


(0.07%) 




0.485 = 


t 0.011 


(2.27%) 


l0g(Me„v/M.) 


-4.39 = 


t 0.10 


(2.23%) 


R/Re {M„ g) 


0.147 = 


t 0.002 


(1.57%) 


L/L^ (r,ff, R) 


28.6 = 


t 1.7 


(5.79%) 


Mv (g, Teff, M.) 


4.44 = 


t 0.05 


(1.13%) 


d (V, Mv) (pc) 


1461 = 


t 66 


(4.52%) 


^ From spectroscopy 



magnitude larger than the measurement uncertainty on the peri- 
ods. This is attributed to inadequacies in the 2G models, which 
can undoubtedly be improved upon. Work on th e so-called third- 
generation models dBrassard & Fontainel2008l) is on-going, and 
should help address at least some of the outstanding issues. 

The solution obtained from asteroseismology leads to a nat- 
ural determination of the three variable input parameters log g, 
M, and log q{li), while the effective temperature is known from 
spectroscopy. On the basis of these quantities we can derive a 
number of secondary parameters: the stellar radius R (as a func- 
tion of g and M»), the luminosity L (as a function of T^ff and R), 
the absolute magnitude My (as a function of g, Teff and M» in 
conjunction with detailed model atmospheres) and the distance 
from Earth d (as a function of apparent magnitude and My', we 
assumed y= 15.26+0.05). As detailed in Section 2, the rotation 
period can tentatively be set to > 2.63 d from the rotational 
splitting observed. The derived parameters are listed in Table |5] 
together with the inter nal 1-cr e rrors computed according to the 
recipe given in Charpi net et al.l (2005a) for the three input pa- 
rameters, and propagated through for the secondary quantities. 
All the errors are purely statistical, and almost certainly under- 
estimate the true uncertainties arising from systematic effects. 



4. Conclusion 

In this paper we presented the asteroseismological analysis for 
a tenth rapidly pulsating subdwarf B star, the relatively recently 
discovered EC 09582-1 137. On the basis of ~ 30 hours of SUSI2 
time-series photometry we uncovered 5 independent harmonic 
oscillations as well as two periodicities interpreted as the rota- 
tionally split components of a frequency multiplet. The first har- 
monic of the dominant oscillation was additionally uncovered, 
albeit at an amplitude below the imposed detection threshold. 
We also obtained a high S/N low resolution spectrum in order 
to infer the target's atmospheric parameters with some accuracy. 
Using the observed oscillations as input, we conducted an astero- 
seismic search in parameter space to find the model that could 
optimally account for the period spectrum of EC 09582-1137 . 
Given that we had only five observed periodicities as input we 
placed some constraints on the mode identification, requiring the 
two dominant pulsations to correspond to radial or dipole modes. 
It was then possible to isolate a well-defined optimal model that 
could reproduce all the observed oscillations simultaneously to 
within 0.6%. As has become standard in the asteroseismology 
of EC 14026 stars, the solution is in accordance with the spec- 
troscopic estimates of the atmospheric parameters, and all the 
observed modes are predicted to be excited from non-adiabatic 
calculations. 



Judging by the optimal model identified, EC 09582-1 137 ap- 
pears to be a very typical member of the EC 14026 pulsator class 
in terms of fundamental parameters. The values of log ^ and Teff 
place it in the middle of the instability strip, and its mass is close 
to the canonically expected value. Our target shows no spec- 
troscopic indication of a companion, and appears to be rotating 
slowly, with a rotation period of the order of a few days. Given 
that sdB stars in close binary systems with Pt,i„ < 0. 6 days have 
been shown to exhibit binary-s ynchronous rotation jGeier et alj 
120081 : Ivan Grootel et alj|2008 ai). such objects should display ro- 
tational splitting corresponding to the binary period in their fre- 
quency spectra. Despite a sufficient time baseline, we find no 
indication of such a short rotation period in our photometry, and 
conclude that, if EC 09582- 11 37 forms part of a binary system, 
the latter has a period longer than ~ 0.6 days. It is therefore a 
good candidate for a single star or a component of a relatively 
wide binary with an unseen companion, but this remains to be 
clarified from radial velocity measurements. 

The fact that the work presented here constitutes the tenth as- 
teroseismological analysis of a rapidly pulsating subdwarf B star 
has in our opinion demonstrated beyond doubt the basic validity 
of our 2G models and the GA approach. Of course, there is still 
room for improvement on the theoretical side. Shortcomings in 
the cuiTent models are indicated by the relatively poor disper- 
sion of the asteroseismic period fit when compared to the obser- 
vational accuracy of the periodicities. The deficiencies of the 2G 
sdB star models become even more apparent when studying the 
slowly pulsating hot subdwarfs. These oscillators exhibit grav- 
ity modes that probe deeper inside the star than pressure modes, 
and are as such sensitive to the inner layers not modelled ac- 
curately by the 2G models. Work on so-called third generation 
models incorporating nuclear processes is ongoing, and is hoped 
to improve the accuracy of asteroseismic solutions. 

On the observational side, much remains to be done. 
Although we were allocated 5 nights of observations on a 
medium-size telescope using a highly sensitive imager, we de- 
tected only 5 independent harmonic oscillations. As a result of 
this, we were forced to make some assumptions as to the mode 
identification based purely on visibility arguments. A longer 
time baseline allowing us to analyse the rotational signature 
and/or multi-colour information would yield partial mode iden- 
tification, making such assumptions unnecessary. Moreover, an 
increased sensitivity would provide additional secure frequen- 
cies, strengthening the robustness of the asteroseismic solution. 
Possible avenues for future observational campaigns include 
space telescopes designed for asteroseismology such as Kepler 
or Corot, and networks of small telscopes providing world-wide 
coverage. An alternative is to obtain shorter high quality data sets 
on large telescopes, however there is currently little instrumen- 
tation available for fast time-series imaging. We remain hope- 
ful that this will change in the near future, e.g with the second- 
generation suite of instruments on the VLT or the new instru- 
mentation for GTC. 
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Fig. 5. Fourier transform of the entire data set zoomed in to the 
5-9 mHz range (the spectrum outside this range is consistent 
with noise). The lower transforms show the successive steps of 
prewhitening as indicated for frequencies detected above 4 cr. 
They have been shifted arbitrarily along the y-axis for visualisa- 
tion purposes, however the amplitude scale is the same as for the 
original Fourier transform. 
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Fig. 7. Nightly variation in amplitude (top panel) and phase (bot- 
tom panel) of the 143 s complex (dots with error bars). The 
dashed (continuous) curves show the amplitude and phase vari- 
ations as simulated for the beating action of the triplet structure 
/3-, /}, /3+ (doublet structure /s, /3+). The amplitude/phase and 
the associated error derived for the full data set is indicated by 
the horizontal dotted hne. 
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Fig. 10. Comparison of the observed period spectrum of EC 
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riod spectrum of the optimal model (dotted line segments) in the 
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Fig. 9. Expected (solid curves) and observed (dots with error 
bars) correlations between the effective temperature, the frac- 
tional hydrogen envelope mass, and the total mass (thick bars 
projected onto the theoretical mass relations) for the 10 sdB pul- 
sators so far submitted to asteroseismology. EC 09582-1137 is 
indicated by the darker bar. 



